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Forces in Oscillating Liquid . 259 

which can be produced in an ice-crystal at right angles to the optic axis very 
closely obeys the laws of viscous flow. If this view should be sustained, it is 
not only true, as Professor H. A. Miers shows, that a liquid of low viscosity 
may have a crystalline structure, but also that, as in the case of crystalline 
ice, a solid may be liquid along one plane only. 

Prom the equation to the curve on Diagram I, the viscosity of an ice- 
crystal in a direction at right angles to the optic axis is about 2 x 10 10 at the 
freezing point. At this temperature the viscosity of a glacier is about 
12,500 x 10 10 , or 6250 times as great. The effective viscosity of a glacier is 
therefore due in a great measure to some other consideration. The optic axes 
of glacier grains are at all angles and they lock each other. 

I am much indebted to Mr. P. H. Parr for the assistance he has rendered in 
calculating the results from the figures given. 


Vortices in Oscillating Liquid. 

By Lord Rayleigh, O.M., Pres. R.S. 

(Received July 25, 1908.) 


In a paper “ on the Circulation of Air observed in Kundt/s Tubes, and on some 
Allied Acoustical Problems,”* I applied the equations of viscous incompressible 
fluid to show that the effect of the bottom of the containing vessel was to 
generate permanent vortices in the vibrating fluid. It was remarkable that 
the intensity of the vortical motion, when fully established, proved to be 
independent of the magnitude of the viscosity, so that the effects could not 
be eliminated by merely supposing the viscosity to become extremely small. 
The expression found for the vortices was simple. The horizontal component 
u of the primary motion near the bottom being u = u 0 cos kx cos nt, the 
component velocities of the vortical motion are 


u = 


, Uq 2 sin 2 kx 


-2ky 


(1-2 ley), 


, u 0 2 cos 2 kx 

v— 4 


~ 2k y 2ky, 


y being measured upwards from the bottom, and V ( = n/k) the velocity of 
propagation of waves of the length in question. According to these 
* ‘Phil. Trans.,’ vol. 175, p. 1, 1883 ; ‘Scientific Papers,’ vol. 2, p. 239. 
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260 Lord Rayleigh. [July 25, 

expressions, the vortical motion is downwards over the places where u has its 
greatest alternating values. In the case of water contained in a tank and 
vibrating in its simplest mode, the theoretical motion is downwards in the 
middle and upwards at the ends. To guard against misinterpretation, it may 
be well to add that quite close to the bottom the motion, as calculated, is of 
a quite different character. 

In a recent paper,* Mrs. Ayrton has examined, with much experimental 
skill, the vortices arising when water oscillates in a narrow tank, and has 
obtained results which differ somewhat widely from what are indicated in the 
above formulae. Near the bottom, and especially when the depth is small, 
there are indeed vortices of this character; but, in general, the most conspicuous 
feature consists of vortices revolving in the opposite direction, the water 
rising in the middle of the tank and falling at the ends. The first thought 
that occurred to me was that Mrs. Ayrton’s vortices might be due to defect of 
freedom in the surface, such as might be supposed to arise from a greasy 
film opposing extensions and contractions ; but in some experiments that I 
tried, the vortical motion did not seem to be much influenced by cleansing 
the surface, and the question was suggested as to whether the free surface 
itself might not originate vortices in somewhat the same way as the bottom 
does, and more potently on account of the greater velocities of the primary 
motion there prevailing. I do not remember whether I had any clear view 
on this question when I wrote the former paper. Yortices originating other¬ 
wise than at the bottom were ignored, but I may not have intended to exclude 
their possibility. 

The present paper consists mainly of an attempt to answer the question 
thus suggested. The investigation is limited to the case of deep water, and 
even then is rather complicated. If the calculations are correct, we are to 
conclude that a free surface does not generate vortices of this kind, at least 
if we suppose the viscosity small and include only the square of the motion. 
How then are the conspicuous vortices, observed by Mrs. Ayrton, to be 
explained ? One might attribute them to the ends of the tank, acting much 
as the bottom does in my former investigation. In the latter case the seat of 
the forces is near the parts lying midway between the middle of the tank and 
the ends (Jtx = + ^tt), and their effect is to push the neighbouring fluid in the 
direction away from the place of greatest motion (kx = 0). A like action at 
the ends of the tank would push the neighbouring fluid downwards, and thus 
generate vortices revolving in the observed direction. An objection to this 
view lies in an observation by Mrs. Ayrton on water oscillating in a long 

* “On the Non-Periodic or Residual Motion of Water moving in Stationary Waves,” 
‘Roy. Soc. Proc.,\A, vol. 80, p. 252, 1908. 
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tank with several subdivisions (p. 255). Since no solid walls are situated at 
the intermediate nodes, it may he thought that the action at the ends would 
be insufficient to establish the whole system of vortices. Probably this would 
be so. But another influence acting in the same direction may arise from the 
faces of the somewhat narrow tanks employed. It would seem that at the 
nodes, lex = (2m+1) 7 r, the friction in moving up and down along the 

faces might have the same general effect as the up-and-down motion along the 
solid walls which constitute the ends. 

But I must confess that some observations made with the help of 
Mrs. Sidgwick were not favourable to this view. The ends of the tank 
were eliminated by using the annular space included between two coaxal 
cylinders, A, B (beakers), but the vortical motion did not seem to be 
diminished. The insertion of a strip of glass C held vertically across the 
annulus, and thus virtually restoring one end, did not make much, if any, 



difference. This experiment seems to exclude the explanation depending 
upon the action of the ends , and that which would attribute the effect to the 
faces is difficult to reconcile with the highly localised character of the effect, 
which at first seems to be limited to the immediate neighbourhood of the ends. 
Can it be that the true explanation would require the retention of terms of 
higher order than the square of the motion ? 

I may mention that on more than one occasion I have witnessed the 
reversed movement described by Mrs. Ayrton on p. 259, “ as if a set of water 
springs had been wound up, and now proceeded to unwind themselves.” 
I presume that the spring depends upon gravity, acting upon unequal 
densities in the fluid, due either to temperature or to variations in the 
amount of powder held in suspension. Another possible explanation would 
lie in the effect of surface contamination. 


In the usual notation the equations of motion in two dimensions are 

1 dp du 
p dx ' 

1 dp 

p cly~ I '~2 dy 


dt 


+ z^y%- 


d (u? + v 2 ) (du 
A —- v\ 


dx 


dy 


civ 

dx 


dv 

= ~17+vV 2 v 


1 d (to 2 + v 2 ) /du dv\ 

Xy ~ U \Xy~Txr g ’ 


( 1 ) 

( 2 ) 
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where y is measured vertically upwards, and v(= fi/p) is the kinematic 
viscosity. Since the fluid is supposed to be incompressible, 


du dv 
dx dy — J 

(3) 

or what is equivalent, 


u = dtyfdy, v = — dtyjdx, 

(4) 

^ being the stream-function. 


In virtue of (4), we have in (1), (2) 


l> 

II 

1 

$ 1 

(5) 


It may be well to commence with the comparatively simple question of 
stationary waves, carried to the second order of approximation, when viscosity 
is neglected. If we eliminate^? from (1) and (2), putting at the same time 
v — 0, we find 



n 


cly 2 ^ dy 2 ± 
dx dy 


( 6 ) 


As was to be expected from the general theory of a frictionless fluid, the 
solution of (6) to any order of approximation is 

VV = 0. (7) 

We now assume that the motion is periodic with respect to x —proportional, 
say, to cos hx, so far as the first approximation is concerned. Thus for this 
approximation we take 

'i/r = Ae ky cos hx cos nt, (8) 

the term in e~ ky , otherwise admissible, being excluded by the consideration 
that all motion must vanish when y = — oo, inasmuch as the fluid is 
supposed infinitely deep. Corresponding to (8), 


dtyjdy = hAe ky cos hx cos nt 

(9) 

—d^jdx = hAe ky sin hx cos nt , 

(10) 

u 2 -f- v 2 = h 2 A 2 e 2ky cos 2 nt. 

(11) 


In virtue of (5), (7), the equations of pressure become 

1 dp __ du ^ d (v? -f- v 2 ) 
p dx~ dt 2 dx 

p dy~ ~~ dt~ 2 dy 

and thus 


rpjp = ~gy+nAe ky sin hx sin nt — cos 2 nt +f (t), 


( 12 ) 
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where f(t) denotes a function of t which is arbitrary so far as the differential 
equations of pressure are concerned. The pressure at the surface is to be 
found from ( 12 ) by putting y = 77 , where rj is the elevation of the surface 


at the point in question. The relation between rj and u , v is thus required 


accurately to the second order of small quantities. 

The differential relation* is 

dv di 7 n m 

dt - V ~ U dx' ( ^ 

To the first order we have drjjdt equal simply to the value of v at the 
surface, so that by ( 10 ) 


rj = hn l A sin kx sin nt, 


(14) 


the origin of y being in the undisturbed surface. This value of r\ may be 
used in the small terms of (13), and thus to a second approximation 



(15) 


where F (x) is an arbitrary function of x of the second order of small 
quantities. 

We are now prepared to substitute for y its value y in ( 12 ). In the 
principal term we must use the complete value of 7 ] from (15). In the second 
term, already containing A as a factor, the first approximation for rj suffices, 
while in the third term we may put y = 0 . The third term thus becomes a 
function of t only, and may be regarded as cancelled by/(£). We find 




(16) 


In free vibrations the condition to be satisfied at the surface is p = 0. The 
annulment of the term in nt requires that 


(17) 


and the same well-known relation suffices to annul the term in 2 nt. The 
surface condition is satisfied without any addition to if besides satisfying 
(17) we identify g F (x) with — ^/c 2 A 2 cos 2 k%. Thus, writing A' for kA/n, we 
obtain, as the complete value of 


rj = A' sin hx sin ^ + pA' 2 cos 2kx cos 2nt— pA' 2 cos 2 hx. (18) 


We now proceed with the consideration of the problem when viscosity 
* Cf Lamb’s 4 Hydrodynamics/ § 10. 
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is retained. Eliminating y> from (1) and (2), we get, with use of (3) 
and (5), 

v 4 t _li v u = ! ( W + !^, (19) 

v dt v dx v dy 

The terms on the right hand of (19) are of the second order in the ampli¬ 
tude of vibration, and thus for the first approximation we have simply 




The solution of (20) may be written 

yjr = ^1 + ^ 2 , 


where 


V 2 fi = 0, 


v2 -^)* 2 = 0 - 


( 20 ) 

( 21 ) 

( 22 ) 


We now introduce the suppositions that in the first approximation yfri, ^]r 2 
are proportional to cos kx, and also to e int . The wave-length along x is 2nrjk> 
and the period r is 27 r/n. The equations (20) now become 

(|- F )+‘ = 0 ' w 


tp 


0, 


(23) 


by which yjn and are to be determined as functions of y. If we write 

Id 2 = k 2 +injv, (24) 

we have, as the most general solutions of (23), 

'xjri = AV^ + B'e""^, = C'e*'v + D'e-*'v; (25) 

but if the real part of Id is taken positive, the terms in B' and D' 
are excluded when the fluid is treated as infinitely deep. Thus for our 
purpose 

yfn = A e k y cos kx e int , yjr 2 = G e y v cos kx e in \ (26) 

where A and C are now absolute constants, real or complex. We shall 
presently find it convenient to suppose C real. From (4) we now find 

n = (kA.e k v -f Id Ged' y ) cos kx e int , (27) 

v = ( kAe k v -f kC e h>y ) sin kx e int ; (28) 

and since to this approximation dyjdt = v, 

Tf = (k/in) (A + C) sin kx e in K. (29) 

If we omit the terms of the second order in (1), we find 

1 dp __ d Idyp'i ^ dyfrd ^ 

p dx 

in virtue of (22) ; and in like manner 

1 dp __ d__ dtyj 
p dy dx dt 


dy\ d,' + ~-if) +V Ty^ + V2 * !) = “I 


-9- 










Downloaded from rspa.royalsocietypublishing.org 


1908.] 


Vortices in Oscillating Liquid. 


265 



(30) 


the expression to be integrated being a perfect differential by (22). 
Applying (30) to the present case, we find 

pj p = — gy — inAe ky sin kx e mt . 


(31) 


At the surface we are to suppose y — % and rj may be neglected in the 
term already multiplied by A. Thus at the surface 



(32) 


Now that we have to reckon with viscosity, the stress conditions at the 
surface can no longer be expressed merely by p. In the usual notation, 
applicable also in the theory of elastic solids, p xx , p yy denote normal 
tractions across faces perpendicular to x and y respectively, while p^ denotes 
the tangential traction which acts parallel to y across the face perpendicular 
to x, or the equal traction parallel to x across the face perpendicular to y. 
The expressions for # these tractions are* 


Pxx — —p + 2/A^~, Pyy — —p + 2p, ( ~~, 


(33) 



(34) 


When viscosity was neglected, we were able to suppose that the surface 
of the fluid was entirely free from imposed force. Under such circumstances 
the vibrations of a viscous fluid could not be maintained. If n is to be real, 
some maintaining forces are necessary. We will suppose that these forces 
are exclusively normal in their character, and accordingly make p^ = 0 ; for 
it is to be observed that in the present approximation a direction parallel to 
the surface may be identified with the horizontal. By (27), (28), (34) we find, 
making y = 0, 


2k 2 A + (k 2 + k' 2 )C = Q, 


(35) 


as the condition of no tangential force at the surface. Or, if we substitute for 
k' 2 its value from (24), (35) becomes 



(36) 


For the normal traction at the surface we have from (33), 
n Pyy!P = A {i (n 2 — kg) -f 2 vk 2 n} sin kx e int 
+ C { — ikg -f 2 vkk’n) sin kx e int , 

* Cf. Lamb’s ‘ Hydrodynamics, 5 § 314. 


(37) 
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and by (36) the expression for 7) in (29) may be written 

C 

7} = — —— sin kx e int . 

2kv 


(38) 


These equations constitute the complete symbolical solution of the problem 
of infinitely small stationary waves maintained by purely normal surface 
pressures in a fluid of any degree of viscosity.* 

In preparing to pass to real quantities, it is simplest to suppose 
C = — 2 kv, so that 

= sin kx e? nt . (39) 


A is then given explicitly by (36), and on substitution in (37) we get 

Pyyjp = {g—n 2 lk-\-4:ivnk-\-4zv 2 k 2 (Jc—k')\ sin kxe int . (40) 

The passage to real quantities is now only complicated by the term in 
Id. If the viscosity be small, this term may be omitted, and we may take 

Pyyjp = {g—n 2 /k + 4dvnk} sin kx e int = (g—n 2 lk) 7} + 4 vk . dy jdt , (41) 

giving the normal traction necessary to maintain the waves represented by 
7] = sin kx cos nt. If n 2 = gk , the part of p yy in the same phase as 7f 
disappears, and 

p y y = 4zpik. d^jdt (42) 


simply. It is to be remembered that p yy is a traction . If, as is usual in 

hydrodynamics, we use pressure ( p '), we see that the pressure has its 

maximum value when the surface is falling fastest, as was to be expected. 

The accurate expression for the real part of (40) may, of course, be 
formed. If in (24) we put 

k 2 = P 2 cos 2a, n/v = P 2 sin 2u , (43) 

then Id = P cos a-hi P sin a. (44) 

It is unnecessary to write down the actual form of the real part of (40).. 
In most applications an approximate value of k' suffices. On account of the 
smallness of i >,njv is very large in comparison with P, that is to say, the 
thickness of the stratum through which the tangential motion can be 
propagated in time t is very small relatively to the wave-length X. We 
may, therefore, usually neglect Id in the equation 

P 4= P + P/P, 

and take simply P 2 = njv . 

Again, (sin a—cos c*) 2 — 1 —sin 2a = ^ldv 2 fn 2 } 

so that the difference between cos a and sin a may often be neglected. 


* Cf. Basset’s ‘ Hydrodynamics, 5 § 520 ; Lamb, loc. cit ., § 332. 


(45) 

(46) 

(47) 
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It appears that the terms neglected in (40) when (41) is substituted are of 
the order v 3/2 


In proceeding to a second approximation we have first to calculate the 
terms of the second order forming the right-hand member of (19), using the 
values found in the first approximation. In these 


' 2 yj/' = v 2 ^ 2 = - e h 'y cos kx e in \ 

v at v 


V Y _ v r , _ - _ 

whence, the imaginary part being rejected, 

= sin kx e v cos a • y sin (nt 4 * P sin a. 3 /) 

CftC v 

= — ”“ C0S ^ 6Pc0Sa,2 ' s i I1 (^ + P s i na -2/ + a )* 

Also, in real quantities by (27), (28), (36), 


-cos^ + ~— sin nt 
2k?v 


u = C cos kx k e ky ^ — ( 

+ P # Pcosa • 2 / cos (nt + V sin a . y+a.)^ 


v = /bC sin kx 


- cos nt -f ~~ sin nt 
2 o 


Hence 


(48) 

(49) 


(50) 


+« p eoB a ■ v cos (rti + P sin «. y) . (51) 

2 p (ud^ 2 yfr/dx + vd^ 2 yfr/dy) 
nkC 2 sin kx cos kx e F cos a • v~ 

= — k e k y {sin (2^ -f P sin a. 3 /) + sin (P sin a. 3 /)} 

-f 2 ^ ^ C ° S ^ s * n ^ ’ 2 /) — cos( 2 ?i£-fPsin a. 3 /)} 

4 -P e k v {sin (2nt + ¥ sin a ,y + a) -f sin(P sin a . y-f a)} 


wP ^ 

{cos (P sin a . y-\-a)— cos(2^-f P sin a . 3 / 4 - a)} 

(52) 


" 2k 2 v 
-2P<? PeOB “-2'sin. a. 


By (19) the equation with which we have to deal is 


y 4 \^ —i yV 2 ^ — ■— x gVcosa.y x right-hand member of (52). ( 53 ) 

V Clo ^xV 


It will be observed that in (52) or (53) the terms on the right, regarded 
as functions of t, are either independent of l or circular functions of 2nL 
Corresponding terms, proportional also to sin 2 kx, will appear in the direct 
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268 Lord Rayleigh. [Jnly 25 , 

integral of (53), and in addition we must include a “ complementary 
function ” representing, so far as required, the complete integral of (53) 
when the second member is made equal to zero. This part contains the 
terms of the first approximation (26), which now represent themselves; and 
we must also be prepared to admit terms of the second order proportional 
to sin 2 kx, and either independent of time or involving 2 ni. In the former 
case the differentia] equation reduces to 

v‘+= (!,-«■)+= °. (54) 

giving as the solution applicable to deep water, 

i/r = sin 2 kx (H -f- K y) e 2ky , (55) 

where H and K are constants. A term in K would represent vortices of 
the kind found in the former paper to arise from the action of the bottom 
(when the liquid is not too deep); and one of the principal objects of 
the present investigation is to ascertain whether these terms occur as the 
result of the conditions operative at the free surface. It may be recalled 
that though such vortices could not arise in an ideal frictionless fluid, their 
magnitude, when fully established, may be independent of the amount of the 
friction. In view of the complication of the problem it must suffice to 
limit the investigation to the case of small viscosity, the question being 
whether vortices can be maintained when the viscosity is reduced without 
limit. 

On the right of (52) there are nine terms, of which five are independent 
of t. But they are not of equal importance. Since P is of the order v~*, the 
leading term is the seventh in order. So far as this term is concerned, 






n 2 PC 2 sin 2 kx 
‘ 8k? " 


. e ky 6 Pcos a • y cos (P sin a . y -f a). 


It is now convenient to revert to complex quantities, regarding 
gP cos a. y cos (p sin cL.y + ct) as the real part of e k 'y +ia . Thus 


giving 


' cl 2 a 7 2 \ 2 » &VC 2 sin2fcr 

__4Pj* =- 


k'n 2 C 2 sin 2 kx 


g( h f +k) y' 


(56) 


8^ 3 {(/ S , +*) 2 -4F} 2 " 

To form the surface condition, representing the evanescence of tangential 
force, we shall require the expression of 

dM + dv = dPjr_d^Jr /(^ + M2 + 4 p\ ^ 
dy dx dy 2 dx 2 
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in which, since C 2 is already involved as a factor, we may put y = 0. Thus. 

(57) 


du dv _ kjrd C 2 sin 2kx { (k' + 7g) 2 +47g 2 } 
dy dx 


8kv* {(Id + k) 2 -4:k 2 } 2 

from which the imaginary part is to be rejected. 

In tracing the value of (57) as v diminishes, we get ultimately 

n 2 G 2 sin 2kx 


(58) 


8 kk'v 3 

In this G 2 /v 2 is of the order A 2 , and Id is of order v~*, so that (58) is of 
order v~^A 2 , becoming infinite in comparison with A 2 , as v diminishes without, 
limit. We shall find, however, that this infinite term is compensated by 
another, to be brought forward later. For our purpose we must retain alt 
terms which do not vanish with v in comparison with A 2 Hence, witfci 
sufficient approximation, 

k' {(k'+k) 2 +4k 2 } l _ V „ 1 2k __ 1 2k 

v{(ld + k) 2 -4k 2 } 2 v(k' + k) 2 vk' vk' 2 vld 


= J_ + 2fifc 

in + vk 2 vk' n ’ 


v(k' + k) 2 

of which the second term, being purely imaginary, is to be rejected. Thus 

n 2 G 2 sin 2 kx n 2 G 2 cos a sin 2 kx 


du , dv -| , n 

-— H—-— = real part or 
dy dx 


(59) 


8kk'v 3 8kv 3 P 

Eeferring back to (52) and having regard to the orders of the various terms 
in respect of v, we see that the only other term independent of t which needs 
be retained is the third in order of arrangement. From this term we obtain 
in like manner, 

. _ n 2 G 2 sin 21cx e (*'+*)y 
Y 8v 3 4**> a * 

and when y = 0 , 

du dv _ m 2 C 2 sin 2 kx {Id + k) 2 + 47 c 2 _ n? C 2 sin 2 kx . x 

dy dx 8v 3 {(/c' + /c) 2 - 4 /c 2 } 2 “ 8vW* (bl ' 

with sufficient approximation. 

In (61), vk 12 , so far as it need be retained, is a pure imaginary, and 
accordingly there is no contribution from this source. We are left therefore 
with (59) as the complete contribution of the direct integral of (53) when 
v = 0 , so far as the terms independent of t are concerned. 

In a similar manner we may treat the terms in 2 nt. We find for the. 
leading term 

V 4 f -- 4 V 2 f = real part of «<*+*>,^. 

v dt 8kv* ? 

so that, if we finally reject the imaginary part, 

du dv __ n 2 G 2 sin 2kx e 2int / 1 2ik\ 

dy dx 8kv 2 \ vk! n ) * 


(62) 
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when y = 0. And for the only other term which it is necessary to retain, 


du^dv __ nC 2 
dy dx Sv 2 


sin Zkx sin 2 nt. 


(63) 


We have now sufficiently complete expressions for du/dy -f dv/dx so far as 
it results from the direct integral of (53). But to these we have to add 
terms arising from the complementary function. In this there must, at any 
rate? be included the terms in nt found in the first approximation. From 
(27), (28) we get, when y is small, 

■ — + $ = cos lex. e int { W A + (k 2 + k' 2 ) C + 2k 2 A . ky + (k 2 +k' 2 ) C. k'y}, (64) 
dy dx 

and in this we are to substitute for y the value appropriate to the surface. 
In the first approximation the terms containing y were neglected, and the 
surface condition gave 


2n 2 A.-\-(l<?-{-k' 2 ) C = 0. 


(65) 


This relation must still hold, approximately at any rate. Using it in the 
small terms of (64), we get 


(k' 2 + k?) G (k' — k) y cos kx e int , 


in which we may neglect k in comparison with k'. Passing to real quantities, 
we find 



The real value of y 9 or rj from (38), is 

7 ] = — sin kx cos nt 9 
2kv 


( 66 ) 


,so that this part of (64) becomes 


-j- ^ == sin 2 kx {sin (2nt + a) + sin a}. 

dy dx okv 2 


sin 2 kx {sin (2 nt H- a) + sin a}. 


(67) 


Thus for the terms independent of t we get, from (59), (67), 



( 68 ) 


The factor within braces vanishes with v, as maybe proved from (43), 
the two infinite terms cancelling without finite residue. Thus altogether 
dujdy + dvjdx vanishes with v so far as the terms independent of t are 
concerned. 

The terms in 2 nt are of less interest. We find, in the same Way, when 
v is diminished without limit for the complete value of duj dy + dvj dx at the 
surface, 


n C 2 sin 2 kx sin 2 nt 

_ ___ 


( 69 ) 
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But we are not yet in a position to apply the condition which must be 
satisfied at the surface, viz., that the tangential force shall there vanish; for 
we must remember that the surface can no longer be treated as parallel to 
y = 0. If 6 be the angle which the surface at the point under consideration 
makes with y = 0, the formulae, of transformation are 


Px'x' = COS 2 0pxx + sin 2 6 p yy + sin 2 6 p xy , 
p y > y > = sin 2 6 p xx + cos 2 6 p yy — sin 26 p xy> 

Px'y> = (cos 2 6 —sin 2 6) p xy + sin 6 cos 6 (p yy —p xx ). 

For the present purpose 6 may be regarded as a small quantity, equal to 
; ±dr)ldx , whose square may be neglected, and we may take 

Px'y' = Pxy + 6 ( p yy p x x)> (70) 

It is p x >y> s and not p xy , which is to be made to vanish to the second order of 
the vibration. 

The expressions for p XX} etc., have already been given in (33), (34). Sub¬ 
stituting the values of u,v of the first approximation, we find, whenz/ is small, 


and 


9 IdU—dv' 
^\dx dy 

drj _ 
dx 


2nC . 7 . 

-- sin Jcx sin 


~ cos kx cos nt. 
2v 


nt , 


(71) 

(72) 


It appears, then, that so far as the terms independent of t are concerned, 
there is no difference between p x > y > and p^ and since we have already seen 
that p^ vanishes, it follows that the surface condition of no tangential force 
is satisfied to a second approximation, without the addition of any further 
terms (55), such as would represent permanent vortices. Accordingly, no such 
vortices exist. 

As regards terms in 2 nt, we find in addition to (69) another term of the 
same form derived from (71) and (72). In a solution complete to the second 
order these terms would need to be compensated by the introduction of new 
second-order terms in i/r of the form 

^ = (Le 2k v -f M^) sin 2 kx e 2int , (73) 

where k" 2 = 4& 2 + 2in/ v ; (74) 

hut it is scarcely necessary for our purpose to define them further. Neither 
does it seem worth while to express at length the equation of pressure when 
the second-order terms are included. The particular case of no viscosity 
already considered illustrates the procedure. Terms in the expression of the 
pressure which are independent of t are balanced by corresponding terms in tj 
not affecting the velocities. 










